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We measure the transverse entanglement of photon pairs on their propagation from the near to the 
far field of spontaneous parametric down-conversion. The Fedorov ratio, depending on the widths 
of conditional and unconditional intensity measurements, is shown to be only able to characterise 
entanglement in the near and far zones of the source. Therefore we also apply a different technique. 
Based on the interference of the two-photon amplitude, our measurement provides direct experi- 
mental access to the Schmidt number of the state. Our measurements show that the interferometric 
scheme always gives the full transverse entanglement of the quantum state. 
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Entanglement is an exciting phenomenon and a funda- 
mental resource in quantum information and quantum 
computation. One convenient source of entanglement 
are photon pairs (biphotons) obtained by spontaneous 
parametric downconversion (SPDC). These photons can 
be entangled not only in discrete variables like polarisa- 
tion or photon number but also position and momentum, 
which are continuous PJH]. The case of continuous vari- 
ables is especially appealing for quantum informational 
tasks because it allows access to a larger Hilbert space 

0. 

Entanglement in the transverse wavevectors of SPDC bi- 
photons can be understood in terms of the famous EPR 
scenario [6] where entanglement is identified by violation 
of the inequality 



(Ax) 2 (Ap) 2 > 
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where (Ax) 2 and (Ap) 2 are the variances of a quantum 
system in position and momentum respectively. The vi- 
olation of Eq. has been measured [8? ] and even 
though it does quantify the amount of transverse entan- 
glement in the biphoton state, measurement in both near- 
position) and far field (momentum) of the source are ne- 
cessary to obtain a value. Closely related to the violation 
of Equation is the Fedorov ratio [9] which is especially 
appealing because it is an entanglement quantifier that 
can be directly measured. Unfortunately the Fedorov ra- 
tio does vary while the state propagates from the near 
to the far field region [10] and even turns to unity at 
some point, indicating no entanglement. Thus the Fe- 
dorov ratio cannot be considered a measure for the full 
entanglement at any arbitrary position. In this work, we 
demonstrate a measurement of the full evolution of the 
Fedorov ratio between those two regimes. Additionally 



we implement a different scheme which has been pro- 
posed [10] to fully quantify the transverse entanglement 
of the biphoton state. This measurement allows direct 
access to the Schmidt number [UHI], a quantity usually 
unattainable in the laboratory but always giving the full 
amount of entanglement. 

Consider the quantum state of a photon pair generated 
by SPDC at the distance z from the centre of the crystal 
as [T2] 
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where p and q are the transverse wavevectors of the sig- 
nal and idler photons and a J (p) and d\ (q) their respective 
creation operators. The properties of the state are gov- 
erned by the two photon amplitude (p, q, z). The most 
widely used operational measure for the entanglement in 
such a system is the Fedorov ratio, which is given by 



R = 
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Here Ap is the standard deviation of the marginal 
distribution of the two photon amplitude P(p,z) = 
f dq |<£ (p, <f, z)| 2 , which we will from here on refer to as 
the unconditional distribution. The width 5p is given by 
the standard deviation of the conditional probability dis- 
tribution P(p\q,z) = |<£ (p, g, z) | 2 at a fixed value of q. 
Analogous expressions hold for p. Another important fea- 
ture of the two photon amplitude is that both transverse 
dimensions are independent so that the two photon amp- 
litude factorises: $ (p, q) = (p xi q x ) (p yi q y ). This al- 
lows us to study the behaviour of the system utilising 
only the spatial degree of freedom in x-direction. 
An intuitive approach to the Fedorov ratio is to think of 
entanglement as strong correlations between very uncer- 
tain occurrences. Thus it is clear that entanglement takes 
a high value for systems where the conditional distribu- 
tion is very narrow (strong correlations) while the uncon- 
ditional distribution is very broad (high uncertainty). 
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The experimental setup is depicted on the right hand side 
of Fig. 0. In order to generate biphotons we focus a cw 
laser with a measured beamwaist of of 245/xm onto a 2mm 
BBO crystal. By means of an interference filter of 8nm 
bandwidth and and apertures positioned at the the cor- 
rect crystal angle, we select the degenerate non-collinear 
type I downconversion process: 404nm — >• 808nm. To se- 
lect a certain part of the angular spectrum, we use slits 
of 30/im width and several mm height. The height of the 
slit leads to an integration over the y direction in both 
conditional and unconditional distributions which does 
not change the ratio between them. Diffraction from the 
slits is compensated by cylindrical lenses before the light 
is fibrecoupled into avalanche photodiodes. For the meas- 
urement one of the slits remains fixed at the maximum 
intensity of the angular spectrum while the other one is 
scanned along transverse direction (in Figure [I] x, per- 
pendicular to the propagation direction z) by means of 
a linear translation stage. We utilise a coincidence elec- 
tronic circuit to record the coincidence-rates as well as 
single count rates of both detectors as a function of the 
transverse displacement. A pair of lenses (f = 500mm) 
is employed to obtain the far field distribution of the 
angular photon spectra (which corresponds to the trans- 
verse momenta) in the focal plane. We apply gaussian 
fits to both the conditional and the unconditional dis- 
tributions and use their respective variances to calculate 
the Fedorov ratio via Equation (J3|. A series of similar 
measurements is performed at various distances z from 
the source up to the image plane of the lenses, where the 
near field (where the coordinate is equivalent to the posi- 
tion quadrature) of the photon distribution is found. The 
three measurements at positions of particular interest are 
shown in Fig. [2] In the same figure one also finds the 
numerical simulations of the two-photon amplitudes in 
the position representation. It is evident that there are 
correlations in the near-field zone and anti-correlations in 
the far field zone while no correlations are observed at a 
certain distance from the crystal. The dependence of the 
Fedorov ratio on the distance is plotted in Fig. [3] Note, 
that because of the effect of the lenses, the far field is ob- 
tained at a smaller distance (focal length) than the near 
field (image plane). It can be clearly seen how the Fe- 
dorov ratio varies as the state propagates through space. 
In the far field we obtain a Fedorov ratio R = 4.0 ± 0.3. 
As we move our detectors towards the near field, the Fe- 
dorov ratio decreases. We even observe a drop to R = 1 
[TU] at a certain distance from the crystal, which would 
indicate no entanglement in this position before the com- 
plete entanglement emerges again in the near field. We 
additionally compare our measurement results with a nu- 
merical simulation (Fig. [3| and find the curve in good 
agreement with our measured data. Since we observe 
that the amount of entanglement increases again after 
the drop and even is completely restored in the near field, 
the entanglement cannot be lost due to decoherence. So 
the question what happened to the 'lost' entanglement 
remains. The explanation, given by Chan et. al. [TO] . 




Figure 1. (Colour online) Experimental setup to determine 
the Fedorov ratio. Measurements are performed at different z 
positions by displacing the detectors. Inset on the left shows 
the modified Mach-Zehnder Interferometer we introduced in 
the signal arm, for the Schmidt number measurement. 



is that upon propagation, the entanglement of the state 
migrates from the wavefunct ion's modulus to its phase 
and back. Thus, since in the intermediate zone between 
near- and far field, the entanglement (or at least parts of 
it) resides in the phase of the state, it is inaccessible to in- 
tensity measurements such as the one we use to determine 
the Fedorov ratio. This prediction is clearly confirmed by 
our measurements. As shown by Tasca et. al. p~3j[14], 
application of different fractional fourier transforms in 
signal and idler channels allows one to reveal intensity 
correlations in the intermediate zone. This is similar to 
the generalised quadrature measurements performed in 
[I5j [16] and requires specific lens configurations for every 
distance z in both signal and idler channels. In the fol- 
lowing we present a more general method, which works 
at arbitrary distances with the same configuration. 
In this method, the measured quantity is the Schmidt 
number K, which is connected to the effective num- 
ber of modes. Any pure state of a composite sys- 
tem AB can be expressed in a certain basis so that 
= Yli VK\^a)\^b)- The number of non-zero ele- 
ments Xi required to express the state vector in terms of 
the Schmidt basis \ia)\^b) is directly related to the lack 
of separability of the state. For infinitely large Hilbert 
spaces the Schmidt number is defined as K = 1/ ^ (A?) 
[T7] . Unlike the Fedorov ratio, the Schmidt number will 
always give the full amount of entanglement of a system 
by definition. The drawback however is that K is usually 
not considered an operational measure, meaning it can- 
not be easily obtained by a direct measurement. Never- 
theless we were able to implement a setup which allows 
us experimental access to the Schmidt number. Since 
the entanglement can manifest itself also in the phase 
of the state, an interferometric setup was required. Un- 
like other schemes, our setup enables one to measure the 
full amount of entanglement without the need to perform 
measurements at specific positions, such as the near and 
far field [18]. Our scheme allows us to obtain the com- 
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Figure 2. (Colour online) Measured conditional (red) and un- 
conditional (blue) photo n distributions with simulated two- 
photon amplitudes, 2(a) in the far field, 2(b) in the interme- 
diate field where R = 1, and 2(c) in the near field. 
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Figure 3. (Colour online) Measurement of the Fedorov ratio 
from far- to near field: Points show measurements, the solid 
line is a numerical simulation. 



plete entanglement from just a single measurement in any 
arbitrary position. Following a proposal by Chan et. al. 
[TO] , we introduced a Mach-Zehnder like interferometer 
in the signal channel. A sketch can be found in the inset 
on the left of Fig. [T]The interferometer contains one dove 
prism in each arm, one of which is rotated by tt/2 with 
respect to the other. Due to internal reflection within 



the prisms the beam in one arm of the interferometer 
is 'spatially mirrored' in one direction while the beam in 
the other arm is 'mirrored' orthogonally compared to the 
first. As this is equivalent to the inversion operation, one 
can say the beam is overlapped with its inverted copy 
on the final beamsplitter. For a single mode state, the 
whole crossection of the beam is able to interfere with it- 
self and we would theoretically expect 100% of visibility 
between the signals at the outputs dj and b\. With an 
increasing number of modes (which can be understood as 
an increase in the degree of entanglement) however, this 
visibility will drop as different modes will not interfere. 
A more rigorous mathematical treatment [10] shows that 
the Schmidt number is given by 



K = 



(P+ + P-) 



{P+-P-Y 

where P + and P_ are the conditional count rates 



dx 8 dxiP^ (x 8 ,Xi) 
dx s dxiP?i -t (x s ,Xi) . 



(4) 



(5) 
(6) 



Here -P a t a t (x 3 , Xj) is the probability to observe the signal 

photon at the output dj of the interferometer at position 
x s and the idler photon at position X{ in the idler mode 
d|. Analogously Pp%\ (x s ,Xi) describes the joint detec- 

s i 

tion probability between the other interferometer output 
b\ and idler d\ at (x 3 ,Xi). 

Figure [4] shows a numerical simulation of the two prob- 
ability distributions at both interferometer outputs. Ac- 
cording to Eq. Q, the Schmidt number is given by the 
inverse visibility of those coincidence probabilities integ- 
rated over x s and X{. Instead of using two detectors as 
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Figure 4. (Colour online) Simulation of the Two phot on in- 
tensities at both output s of t he interferometer, 4(a) shows 



constructive interference, 4(b) depicts the complementary de- 
structive interference between the interferometer arms. 

indicated in Fig. [T]we measured the visibility in coincid- 
ences by scanning the phase of the interferometer with 
the help of a piezo attached to one of the mirrors. This 
time we position two 200 jam slits oriented along the x dir- 
ection in both signal and idler arms. This is equivalent 
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to an integration along x s and x s in one single meas- 
urement. Hence we are able to determine the Schmidt 
number simply by moving the piezo from the position of 
maximum count rate to minimum count rate in one out- 
put of the interferometer and record the visibility. Both 
the results of the measurements as well as the numer- 
ical predictions are depicted in Fig. [5] Please note that 
the numerical model does not perform the actual Schmidt 
decomposition of the quantum state but rather simulates 
the effect of our measurement method on the state. It 
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Figure 5. (Colour online) Circles are measurements of the 
Schmidt number, the solid red line is the result of a numerical 
simulation of the Schmidt number measurement. For compar- 
ison the numerical simulation of the Fedorov ratio is included 
as well (solid blue line). 

is apparent that the measurement points exclusively lie 
above the numerical prediction, which is not surprising 
because the entanglement we are trying to determine ac- 



tually manifests itself in the lack of visibility. Accord- 
ingly any experimental imperfection always will lead to 
an overestimation of entanglement and consequently a 
slightly increased Schmidt number. 

As mentioned before, the reason for the breakdown of R 
is that upon propagation in free space the two-photon 
amplitude acquires a quadratic phase 

$ (q, z) = ® (q, 0) exp (ikz) exp (-*fj^) > ( 7 ) 

where q is again the transverse component of the 
wavevector k. Henceforth, this leads to a shift of entan- 
glement from the modulus to the phase of the quantum 
state at a certain position. The quadratic phase in Eq. 
^ is mathematically identical to the one acquired by the 
two-photon amplitude due to the group velocity disper- 
sion in an optical fibre [19]. Thus it should be possible 
to obtain similar results in the time-frequency entangle- 
ment of spontaneous parametric downconversion photons 
using a fibre as dispersive medium. Such effects could be 
an issue in QKD experiments where two-photon states 
are sent through very long fibres. 

In conclusion we observed the migration of entangle- 
ment in the transverse momentum of the quantum state 
emitted by spontaneous parametric downconversion from 
modulus to phase as the photons propagate through free 
space. We present both numerical and experimental 
results showing the seemingly vanishing entanglement 
between near and far field of the crystal as given by the 
Fedorov ratio. Furthermore we implement a direct meas- 
urement of the Schmidt number, which unambiguously 
gives the full entanglement of the quantum state. Fi- 
nally we propose an analogue experiment which could 
reveal similar behaviour in the time- frequency picture. 
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